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We use configuration interaction technique to calculate parity non-conservation (PNC) in 
metastable Xe and Hg [proposal of the experiment in L. Bougas et al, Phys. Rev. Lett. 108, 
210801 (2012)]. Both, nuclear spin-independent and nuclear spin-dependent (dominated by the nu- 
clear anapole moment) parts of the amplitude are considered. The amplitudes are strongly enhanced 
by proximity of the states of opposite parity. 
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I. INTRODUCTION 

The study of the parity non-conservation in atoms is 
a low-energy, relatively inexpensive alternative to high- 
energy search for new physics beyond the standard model 
(see, e.g. [J-Q). For example, parity non-conservation in 
cesium is currently the best low-energy test of the elec- 
troweak theory J Q • It is due to high accuracy of the 
measurements [4| and its interpretation Q (see also @) . 
Since the cesium result is unlikely to be significantly im- 
proved, the focus of the PNC study in atoms has shifted 
to two important areas: (i) the PNC measurements for 
a chain of isotopes and (ii) the measurements of nuclear 
anapole moment (see, e.g. Q). Most of current or planed 
PNC experiments in atoms consider both possibilities. 

The experiments are in progress at Berkeley for Dy and 
Yb atoms 0, H], at TRIUMF for Rb and Fr atoms [§, 
[l0|, and at Groningen (KVI) for Ra+ ions^lj,]. There 
is an an interesting recent suggestion to measure PNC 
in metastable Xe and Hg In the present paper we 
support this proposal by the atomic calculations. 

The advantages of using Xe and Hg for the measure- 
ments of PNC in a chain of isotopes and the measure- 
ments of anapole moments include: (i) large number of 
stable isotopes of both atoms (maximal difference in the 
number of neutrons is AiV = 12 for Xe and A7V = 8 for 
Hg); (ii) presence of two stable isotopes for each of the 
atoms with non-zero nuclear spin (^^^Xe, / = 1/2; ^'^^Xe, 
/ = 3/2; i^^Hg, / = 1/2; ^oiRg, / = 3/2), these isotopes 
are suitable for the anapole moment measurements; (iii) 
The PNC amplitudes in the Xe and Hg atoms are en- 
hanced due to the high nuclear charge and strong mixing 
with close states of opposite parity. 

An extra advantage comes from the fact that xenon 
and mercury nuclei with non-zero spin have valence neu- 
tron, therefore the nuclear anapole measurements will 
provide the strength constant for the neutron-nucleus 
PNC potential. The anapole moment so far was mea- 
sured only for ^^■^Cs [1] which has valence proton. The 
data for the xenon and mercury would be complimentary 
to those obtained for cesium. 

In the present work we use the configuration inter- 
action (CI) technique to calculate the nuclear spin- 



independent PNC amplitudes caused by the weak nuclear 
charge and the nuclear spin-dependent PNC amplitudes 
dominated by the nuclear anapole moments. The result 
is presented in a convenient form as a sum of two con- 
tributions for different hyperfine structure components. 
This would allow to extract both the value and the sign 
of the anapole moments by comparing the measured am- 
plitudes with the calculated ones. 



II. GENERAL FORMALISM 

The Hamiltonian describing parity-nonconserving 
electron-nuclear interaction can be written as a sum of 
the nuclear-spin-independent (SI) and the nuclear-spin- 
dependent (SD) parts (we use atomic units: h — \e\ ~ 
me = 1): 
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where Gp ~ 2.2225 x 10^^^ a.u. is the Fermi constant 
of the weak interaction, Qw is the nuclear weak charge, 
cr ' 
cr 

nuclear spin, and p{r) is the nuclear density normalized 
to 1. The strength of the spin-dependent PNC inter- 
action is proportional to the dimensionless constant >c 
which is to be found from the measurements. There 
are three major contributions to >f arising from (i) elec- 
tromagnetic interaction of atomic electrons with the nu- 
clear anapole moment [l3| . (ii) electron- nucleus spin- 
dependent weak interaction [14| . and (iii) combined ef- 
fect of spin-independent weak interaction and magnetic 
hyperfine interaction [l^ (see also [l|). In this work we 
do not distinguish between different contributions to k 
and present the results in terms of total x which is the 
sum of all possible contributions. 

Within the standard model the weak nuclear charge 
Qw is given by ISj 

Qw ~ -0.9877iV + 0.0716Z. (2) 

Here A'^ is the number of neutrons, Z is the number of 
protons. 
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The PNC amplitude of an electric dipole transition 
between states of the same parity \i) and |/) is equal to: 
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where d = —e J^i '"'i is the electric dipole operator. To 
extract from the measurements the parameter of the 
nuclear spin-dependent weak interaction k one needs 
to consider PNC amplitudes between specific hyper- 
fine structure components of the initial and final states. 
There amplitudes can be expressed as 
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Here F = J + I, I is nuclear spin. Detailed expressions 
for the reduced matrix elements of the SI and SD PNC 
amplitudes can be found e.g. in Refs. 17 1 and [l^. For 
the SI amplitude we have 



{Jf,Ff\\dsiU,F,) = 

{2Ff + l)i2F, + l){ -^/^ \ 
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= c{Ff,Jf,F,,J,)E'f^. 

Here c{Ff,Jf,Fi,Ji) is the angular coefficient and the 
sum over n, E'j^ does not depend on Ff or Fi: 



^' = E 



(J/||d||n, J„)(7i, Jn\\Hsi\\Ji) 



Ei — En 

{Jf\\Hsi\\n,Jn){n,Jn\\d\\J,] 

Ef-En 

For the SD PNC amplitude we have 
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TABLE I: Energy levels (in cm^^) and g-factors of low states 
of mercury. States considered for the PNC transitions are 
marked as A and B. 

Configu- Calculations Experiment 

ration State Energy g Energy g 



6s" 
6s6p 



6s7s 



Ai 
Aa 

A3 

B 



3r>o 
Spo 
3r>o 
Ipo 

^Si 

^s„ 





38202 
39955 
44812 
53584 
61879 
63399 



0.000 
0.000 
1.480 
1.500 
1.020 
2.000 
0.000 





37645 
39412 
44043 
54069 
62350 
63928 



0.000 
0.000 
1.486 
1.501 
1.015 
2.003 
0.000 



TABLE II: Energy levels (in cm~^) and (?-factors of low 
states of xenon. States considered for the PNC transition 
are marked as A and B. 



Configu- 



Calculations 



Experiment 



ration 




J 


Energy 


9 


Energy 


g 


5^" 










0.000 





0.000 


5p^6s 


A 


2 


64680 


1.500 


67068 


1.501 






1 


66089 


1.214 


68046 


1.206 









74844 


0.000 


76197 


0.000 




B 


1 


75999 


1.314 


77186 


1.321 






1 


76083 


1.853 


77270 


1.852 






1 


77914 


1.024 


78957 


1.022 



The PNC amplitude between different hfs components of 
the initial and final states can be presented in a form 

-EpNc(i^i,i^2) = c{FuF2)E'Q^ [1 + RiFi , F2) , (8) 

where c is an angular coefficient and R is the ratio of 
the spin-dependent to the spin-independent PNC ampli- 
tudes. 

If at least two PNC amplitudes Ei and E2 are mea- 
sured then the value of x can be expressed via the mea- 
sured ratio of the amplitudes Ei / E2 and the calculated 
ratios R of the SD and SI PNC amplitudes. 



El =ciE'Qw{l + Rix), 

E2^C2E'Qw{l+R2>c), 
_ Ci/C2 - E1/E2 

R2E1/E2- R1C1/C2' 



(9) 
(10) 



The ratios Ri , R2 are much less sensitive to numerical 
uncertainties than each of the SD and SI PNC ampli- 
tudes [11]. 



III. CALCULATIONS 

Calculations for xenon and mercury were performed 
with the use of the configuration interaction (CI) method. 
We treat Hg as an atom with two valence electrons and 
Xe as an atom with six valence electrons. Calculations 
for Hg are very similar to what we did before for Hg 
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and Yb [21j. We use the V approximation in which 
initial Hartree-Fock procedure is done for the Hg^+ ion. 
The complete set of single-electron orbitals is constructed 
using the B-spline technique 22|. Core- valence correla- 
tions are included by adding the second-order correlation 
potential S to the CI Hamiltonian in the framework of 
the CI-I-MBPT method [2^. Accuracy for the energies 
is further improved by rescaling the core- valence correla- 
tion operator E (see [13, HII for details). The rescaling 
coefficients are = 0.82 for s-states and Ap — 0.9 for 
p-states. The calculated energies and g-factors of mer- 
cury are presented in TableUtogether with corresponding 
experimental numbers. The g-factors are useful for the 
identification of states. The accuracy for the calculated 
energies is within 1% for majority of the states. It is not 
perfect in spite of the fitting because with use only two 
fitting parameters for all states. 

Similar approach for xenon is problematic due to the 
larger number of valence electrons. We treat all 5p elec- 
trons as valence ones, so that total number of valence 
electrons is six. Using the same technique as for Hg would 
lead to a very large CI matrix. It was suggested in [ll] 
that the hole-particle formalism can be used for the cal- 
culation of the electron structure of xenon. In this case 
only two active particles enter the CI calculations and 
the calculations might not be more complicated than for 
mercury. This assumes that only single excitations are 
allowed from the 5p subshell. However, in our experience 
double excitations are also important. Inclusion of dou- 
ble excitations within the hole-particle formalism hugely 
complicates the problem. Therefore, we use simpler ap- 
proach. We use standard CI technique for six valence 
electrons. Initial Hartree-Fock procedure is done for neu- 
tral xenon (the approximation). Single-electron basis 
states above the core are calculated as the Hartree-Fock 
states in the V'^^^ potential of the frozen core. Many- 
electron basis states for the CI calculations are formed 
by allowing single and double excitations from the 5p 
subshell to the states above the core. 

Accurate treatment of the core- valence correlations for 
xenon within the CI+MBPT method is problematic due 
to large contribution of the subtraction diagrams ^3 . On 
the other hand the core- valence correlations are relatively 
small for xenon compared to valence- valence correlations 
due to large number of valence electrons. Therefore, they 
can be included approximately. 

To simulate the effect of core-valence correlations we 
include into the CI Hamiltonian the core polarization po- 
tential 



6V 



OLc 



(11) 



where ac is the polarizability of the core and a is a cut- 
off parameter. We use a = as and treat ac as a fitting 
parameter. This allows us to fit the energy interval AE — 
84 cm~^ between state B and close state of the same total 
momentum J and opposite parity belonging to the 5p^6p 
configuration. This is important because PNC amplitude 
is very sensitive to this energy interval. We use ac = 



0.554 a.u. for s states of valence electrons and etc = for 
other states. 

The results for energies and g-factors of xenon are pre- 
sented in Table HT] Note that the accuracy is similar to 
what we have for mercury. This is because the the core- 
valence correlations are strong for mercury, they are sig- 
nificantly stronger than the correlations between two va- 
lence electrons. This is a limitation factor for the ac- 
curacy of the calculation. In contrast, the core- valence 
correlations are small for xenon. 



A. Dalgarno-Lewis method for matrix elements 

To calculate PNC amplitude we need to calculate 
matrix elements between many-electron states and per- 
form summation over complete set of many-electron basis 
states (see, e.g. © and ([71). We use random phase ap- 
proximation (RPA) [23, [25I for the matrix elements and 
Dalgarno-Lewis method [26[ for the summation. 

Matrix elements are given by 



£^1™ = {'fyWf + 6Vf\\^, 



(12) 



where SV is the correction to the core potential due to 
the core polarization by an external field /. In present 
calculations / represents either external electric field, SI 
weak interaction or SD PNC interaction. 

Summation over complete set of many-electron states 
is reduced to calculation of the correction 6'^^ to the 
many-electron wave function of the state v due to the 
weak interaction perturbation iJpNC- Then, the PNC 
amplitude is given by 



(13) 



The correction d'i'y is found by solving the system of 
linear inhomogeneous equations 

(il^ff - Ey)6^y = -{Hp^c + (^l^PNc)*.- (14) 

The proposal of [l^ considers the PNC transitions be- 
tween the exited states A and B in Hg and Xe (see Tables 
HI and HI)) . The upper state B in both atoms is very close 
to another state of the same total momentum J but op- 
posite parity. The interval is 8282 cm"'^ for Hg and only 
84 cm~^ in Xe. This is a strong advantage of using these 
transitions from the experimental point of view because 
the proximity of the sates of the same total angular mo- 
mentum but opposite parity leads to the strong enhance- 
ment of the PNC amplitude. On the other hand, this is 
a complication from the theoretical point of view. The 
PNC amplitudes are sensitive to small energy intervals 
where a small error in the calculated energy of the states 
can lead to a large error in the value of the PNC am- 
plitude. To go around this problem we use a stabilizing 
procedure which consists of the following steps. First, 
we use the procedure described above (see Eq. (IT5|) and 
(ITil) ') without modifications. Then we repeat the calcula- 
tions applying the orthogonality conditions for the S'^b 
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to the close state of same J and parity. The contribution 
of the close state is found by comparing the two results. 
Finally, this contribution is added back to the PNC am- 
plitude with the rescaling parameter Ai^thcor/AiJcxp- 

The procedure described above corresponds to the ex- 
act fitting of the energy denominators in ([6]) and (O. 
Therefore, the same results should be obtained if the im- 
portant energy intervals are fitted exactly by rescaling 
the correlation potential S (for Hg) or polarization po- 
tential SV (for Xe) . This is another important test which 
we used in the calculations. 



B. Accuracy of the calculations 

Accuracy of the similar calculations for two-valence- 
electron atoms H g an d Yb were discussed in detail in our 
earlier works [20I l2ll [27j . It was demonstrated that the 
accuracy for transition amplitudes and polarizabilities is 
on the level of 5%. Note that the PNC amplitude is a 
second-order effect similar to polarizability but with one 
electric dipole operator replaced by weak interaction. In 
present work we assume the same 5% uncertainty for the 
PNC effect in mercury. This is supported by the study 
of the limitation factors which are discussed below. 

The uncertainty for xenon is higher due to larger num- 
ber of valence electrons which make it difficult to saturate 
the basis. We assume the 10% uncertainty which come 
from comparing the results obtained with two different 
basis sets, the b-spline basis set and the Hartree-Fock 
basis set. 

The main factor limiting the accuracy of the calcu- 
lations of the PNC amplitudes in xenon and mercury is 
the proximity of the levels of the same total momentum J 
and opposite parity. These states are mixed by weak in- 
teraction and small energy interval between them leads 
to stong enhancement of the PNC effect. This is one 
of the main reason for the choice of atoms and transi- 
tions. However, it represents a challenge for the calcula- 
tions. Even small theoretical error in the energies of the 
mixed states can lead to large error in the PNC ampli- 
tude. There are two ways around this problem. One is 
to fit the energy interval exactly. Another is to isolate 
the resonant contribution from the rest of the amplitude 
and to re-scale it to the correct energy interval. We did 
both things for xenon and found that they give the same 
results. Therefore, for mercury we just fit the energies. 
As long as the problem of small energy intervals is under- 
stood and properly dealt with it didn't contribute much 
to the uncertainty of the results. 

Another important limitation factor for xenon, which 
is harder to deal with, is large number of valence elec- 
trons. It was suggested in Ref. that xenon can be 
treated as a system with one electron and one hole in 
the 5p subshell. This would make it a two-particle sys- 
tem similar to mercury. However, such approach assumes 
that only single excitations from the 5p subshell are con- 
sidered. In our experience full saturation of the basis is 



not possible without inclusion of double and may be even 
triple excitations. This would lead to the CI matrix of the 
extremely large size. Probably, this can be done with the 
use of supercomputers. Another open question for xenon 
is whether and how to include core- valence correlations. 
Incomplete saturation of the basis is the main source of 
the uncertainty of present calculations for xenon. 

Saturation of the basis is not a problem for mercury 
when it is considered as a system with two valence elec- 
trons. However, here we have another problem, strong 
core- valence correlations. The outermost subshell of the 
mercury core, 5d is strongly mixed with the valence 6s 
electrons. This is evident from the presence in the dis- 
crete spectrum of mercury the states with excitations for 
the 5d subshell. The core-valence correlations are in- 
cluded in present work with the use of the second-order 
correlation operator E. Most probably the correlations 
are too strong to be treated accurately within the second- 
order approach. The answer may come from the use of 
some all-order technique similar to what was recently de- 
veloped in Ref. jH. 

Even though some improvement of the accuracy of the 
calculations is possible, there is little chance that it will 
ever match the accuracy for cesium Q. Therefore the 
main focus of the PNC study in Xe and Hg should be 
directed to the measurements of the anapole moments 
and to the study of the ratio of PNC effects in different 
isotopes. 



IV. RESULTS AND DISCUSSION 

The calculated nuclear spin-independent PNC ampli- 
tude for Xe (z-component) is 

£;pNc(A -> B) = 1.76 X IQ-'^" {~Qw/N)ieaB. (15) 

The spin-independent PNC amplitudes for Hg are 

^^PNc(Ai ^ B) = 2.09 X IQ-^" {-Qw/N)ieaB, 
-BpNc(A2^B) = 1.77 X 10-^° {-Qw/N)ieaB, {16) 
-BpNc(A3 ^ B) = 1.25 X 10-^°{-Qw/N)ieaB. 

The difference in the value of the PNC amplitude for 
different isotopes is mostly due to different value of the 
weak nuclear charge Qw Therefore, the amplitudes ([15]) 
and (jl6l) may be used for any isotope. 

Detailed data for both SD and SI PNC amplitudes for 
isotopes with non-zero nuclear spin are presented in Table 
Unifor Xe and Table |TVl for Hg. 



A. Ml amplitudes 



The experimental proposal |12| is aimed to measure 
the PNC optical rotation. The angle of rotation is pro- 
portional to the ratio R — Im{Ep^c/Ml). Therefore, 
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TABLE III: PNC amplitudes (z-components) for the |5p^6s^ ^[3/2]5,Fi) |5p^6s ^[1/2]^, F2) transitions in ^^^Xe and ^^^Xe. 
/ is nuclear spin, F = J + I . E' is given by ((Gj. 



A 


/ 


E' 

units 10~^°jeao 


Fi 


Fa 


PNC 

E'Qw 


amplitude 

10-^°ieao 


129 


1/2 


3.16 


3/2 


1/2 


(l/3)(l + 0.0387>f) 


1.05(1 +Q.0387j<) 








3/2 


3/2 


-l/A/(50)(l + 0.010jf) 


-0.45(1 + O.OlOjf) 








5/2 


3/2 


y(2/25)(l - 0.0226><) 


0.89(1 - 0.0226>f) 


131 


3/2 


3.25 


1/2 


1/2 


-1/^18(1 + 0.0345>f) 


-0.766(1 + 0.0345j<) 








1/2 


3/2 


-1/^90(1 + 0.0252x) 


-0.343(1 + 0.0252j<) 








3/2 


1/2 


1/718(1 + 0.0282x) 


0.766(1 + 0.0282>f) 








3/2 


3/2 


-78/125(1 + 0.0189^) 


-0.822(1 + 0.0189>f) 








3/2 


5/2 


-l/^/375(l + 0.0034>f) 


-0.178(1 + 0.0034>f) 








5/2 


3/2 


77/125(1 + 0.0083>!r) 


0.769(1 + 0.0083>f) 








5/2 


5/2 


-73/70(1 - 0.0072j<) 


-0.673(1 - 0.0072>f) 








7/2 


5/2 


72/35(1 - 0.0220jf) 


0.777(1 - 0.02204x) 



TABLE IV: PNC amplitudes (z-components) 


for the |6s6p 


'P},Fi) - 


-> 6s6p ^Pi,F2) transitions in 


is^Hg and 20iHg. / is 


nuclear spin, F = J + I . E' 


is given by 










A I J 


E' 


F^ 


Fo 


PNC amplitude 










E'Qw 


10""'ieao 


199 1/2 


3.41 


1/2 


1/2 


(1/3)(1 + 0.0084>ir) 


1.14(1 + 0.0084>ir) 






1/2 


3/2 


-7^(1 - 0.0042>f) 


-1.61(1 - 0.0042>f) 


1 


5.31 


1/2 


1/2 


72/27(1 + 0.0538>f) 


1.44(1 + 0.0538>f) 






1/2 


3/2 


71/27(1 - 0.0860>f) 


1.02(1 - 0.0860>f) 






3/2 


1/2 


71/27(1 + 0.0308>f) 


1.02(1 + 0.0308>f) 






3/2 


3/2 


1/76(1 -0.0105>f) 


2.17(1 - 0.0105>f) 


2 


3.71 


3/2 


1/2 


(1/3)(1 + 0.0381>f) 


1.24(1 + 0.0381>f) 






3/2 


3/2 


-1/750(1 + 0.0471>f) 


-0.525(1 + 0.0471>f) 






5/2 


3/2 


72/25(1 - 0.0264>f) 


1.05(1 - 0.0264>f) 


201 3/2 


3.47 


3/2 


1/2 


(l/3)(l + 0.0069>f) 


1.16(1 +0.0069j<) 






3/2 


3/2 


1/75(1 + 0.0028x) 


1.55(1 +0.0028x) 






3/2 


5/2 


-72/15(1 -0.0041x) 


-1.27(1 - 0.0041x) 


1 


5.40 


1/2 


1/2 


-1/754(1 + 0.0171>f) 


-0.735(1 + 0.0171>f) 






1/2 


3/2 


75/54(1 + 0.00362>f) 


1.64(1 + 0.00362x) 






3/2 


1/2 


75/54(1 + 0.0419>f) 


1.64(1 + 0.0419>f) 






3/2 


3/2 


72/75(1 + 0.0607>f) 


0.882(1 + 0.0607>f) 






3/2 


5/2 


(l/5)(l-0.0515>f) 


1.08(1 - 0.0515>f) 






5/2 


3/2 


(l/5)(l + 0.0122>f) 


1.08(1 + 0.0122>f) 






5/2 


5/2 


1/76(1 - 0.0103j<) 


2.20(1 - 0.0103>f) 


2 


3.78 


1/2 


1/2 


-l/v^(l + 0.0385>f) 


-0.890(1 + 0.0385>f) 






1/2 


3/2 


-1/790(1 + 0.0414>f) 


-0.398(1 + 0.0414>f) 






3/2 


1/2 


1/718(1 + 0.0240>f) 


0.890(1 + 0.0240>f) 






3/2 


3/2 


-78/125(1 + 0.0269>f) 


-0.956(1 + 0.0269x) 






3/2 


5/2 


-1/7375(1 + 0.0318>f) 


-0.195(1 + 0.0318>f) 






5/2 


3/2 


77/125(1 + 0.00286>f) 


0.894(1 + 0.00286x) 






5/2 


5/2 


-73/70(1 + 0.00774x) 


-0.782(1 + 0.00774x) 






7/2 


5/2 


72/35(1 - 0.0260>f) 


0.903(1 - 0.0260j<) 



we need to know the values of the magnetic dipole am- 
phtudes for the transitions proposed for PNC measure- 
ments. The most accurate values of the Ml amplitudes 
can be found analytically using experimental values of 
the magnetic g-factors to find the coefficients for config- 
uration mixing. 



This is especially important for the case of mercury 
where numerical calculations of the Ml amplitudes give 
unstable results. The reason for this instability is easy 
to understand. The transitions considered for the PNC 
measurements in mercury are between states of different 
spin (5=1 for states Ai 2,3 and S = for state B). 
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This means that the Ml amphtudes between these states 
vanish in the non-relativistic hmit. In relativistic cal- 
culations the amplitudes are not zero but small. These 
small values are obtained as a result of strong cancelation 
between different contributions. This stong cancelation 
leads to unstable results. 

On the other hand analytical evaluation of the Ml 
amplitudes is simple and produce very accurate results. 
The operator of the magnetic dipole transition (M^ — 
{Lz + 2Sz)fj,B) has no radial part and cannot change a 
principal quantum number in the non-relativistic limit. 
Therefore, the magnetic (^-factors and Ml amplitudes are 
mainly sensitive to the mixing of the states belonging to 
the same configuration. Mixing with other configurations 
normally produces corrections at the 10~^ level [29|. We 
may see this in Table H] and |ll] where factors of " pure" 
sp states with J = 2 and sum of the g-factors for mixed 
states J — 1 differ from the experimental values by less 
than 0.1 %. Therefore, the mixing coefficients for the 
states belonging to the same configuration and Ml am- 
plitudes can be found practically exactly from the known 
values of the 5- factors. Note that we use the calculated 
value of the overlap between the radial wave functions 
Pi/2 and which is close but not equal to 1. For Hg 
it is 0.988. 

We get for the states A and B of xenon 
*A = 1 5^3/265) 

*s = 0.05|5p3/26s) + 0.999|5pi/26s). (17) 



the form 



This leads to the Ml amplitude 



MIab = 1.22^i3 = G.00446eaB. 



(18) 



Using E' from table IHIl we get R = 7.1(7)(35) x 10"** 
for i29xe and R = 7.3(7) x 10"^ for "iXe. Here we 
assume the 10% uncertainty as it has been discussed 
above. These values of Ml and R are close but not 
in perfect agreement to what was found in Ref. p^ : 
Ml = 0.0042eaB, R = 11(3) x 10"^ The reason for 
difference in Ml is not clear. The authors of Ref. [11] 
use slightly different coefficients of configuration mixing 
in (Ull). Their values are 0.062 and 0.998. However, if 
we use these coefficients we get Ml — 0.00444eaB which 
is different from the value All ~ 0.0042eaB presented in 
Ref. 12]. 

The situation is more complicated for mercury. All 
the Ml transitions of interest happen between the states 
with different total spin and vanish in the non-relativistic 
limit since the operator Lz + 2Sz conserves the total spin. 
This leads to the strong suppression. The values of the 
Ml amplitudes presented in Ref. [T2!] are too large for the 
spin-forbidden transitions. No choice of the configuration 
mixing coefficients can reproduce them. 

The wave functions for states Ai^2.3 and B for Hg have 



*Ai = ]6pi/26s) 

- 0.432]6p3/26s) 

*A3 = ]6p3/26s). 

*B = 0.902]6p3/26s) - 



-G.902]6pi/26.s), 
0.432]6pi/26s), 



(19) 



The coefficients 0.902 and 0.432 are chosen to fit the ex- 
perimental 5-factors of the states A2 and B. When pro- 
jections are included all \'^^^^Pjj^) states of the 6s6p 
configuration can be written as 



|Ai 

]A2 

IA3 
IB 



'P?i) 



1 

71 



6s 11 

2 2 



I6s3 



iQp 



2 2 



-0.216]6si i6p3 1 

' 2 2 2 2 



) + 



'P22) 
'Pii) 



0.374]6si_i6p3 3) +0.902]6si i6pi 1), 

' 22 22' ' 22 22' 

]6sii6p3 3), (20) 



-0.451]6si i6n3 1) 

' 2 2 22' 



0.781l6si 



L 6p3 



0.432]6si i6pi 



Here we use (119)) for the expansion. The results for 
Ml amplitudes obtained with the use of these formulas 
are presented in Table |V] together with the values from 
Ref. [ill. Note that there must be Ml <^ p,B for the 
spin-forbidden transitions. This holds for the results of 
present work but not for the results of Ref. [l^- Ta- 
ble |V] also presents the PNC amplitudes and the ratios 
R = Iin{Ep^c/Ml). The numbers include 5% error bars 
according to the estimated uncertainty of the calculations 
which was discussed in section UlIB I Note that the values 
of the Ml amplitudes obtained in the present work are 
practically exact due to the fitting of the experimental 
^-factors. 

The values of the ratios R are larger than in Ref. [T^] 
due to smaller Ml amplitudes. The values of R for mer- 
cury are about an order on magnitude larger than for 
xenon (see above) and about an order of magnitude larger 
than for Tl, Pb and Bi [M^. 



B. Optical rotation 



The angle of optical rotation is given by [33 



0PNC 



47r? 

— (n(a;) 



1)R, 



(21) 



where / is path length in vapor, A is wavelength of laser 
light, uj is its frequency, and n{uj) is refractive index due 
to the absorption. Although the angle is proportional 
to the ratio R — Im{Epjqc/Ml), the small values for 
the Ml amplitudes do not necessarily translate into a 
large angle of rotation. This is because the refractive 
index also depends on Ml amplitude, n{uj) — 1 ^ ]M1]^, 
which leads to </> ~ Ml ■ Ep^c- However, the suppression 
due to the small Ml amplitude can be compensated at 
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TABLE V: Magnitudes of the Ml and PNC amplitudes (reduced matrix elements) and ratios R = /m(i5pNc/Afl) for the 
^Pj — ^ Pi transitions in mercury. Values for Ml from Ref. [12] are also presented for comparison. 



Transition 




Ml amplitudes 




Im{EpKc) 




R 




This work 




Ref. [12] 






is^Hg 


20lHg 




units lO^^eas 




lO^^eas 


fJ-B 


lQ"^°eas 




10-^ 


•"Po" P° 


8.37 


0.229 


14 


0.384 


3.4(2) 3.5(2) 


4.1(2) 


4.2(2) 


^Pi° pr 


7.26 


0.199 


42 


1.15 


5.3(3) 5.4(3) 


7.3(4) 


7.4(4) 


^P2 Pi 


9.94 


0.272 


57 


1.56 


3.7(2) 3.8(2) 


3.7(2) 


3.8(2) 



sufBciently high vapor pressure by the appropriate choice 
of the path length These questions are discussed in 
detail in Ref. [331 ■ Here we just note that the angle of 
rotation per unit length is (j) Ml ■ E'pnc- However, the 
angle of rotation per absorption length (jj ^ R and it is 
large for small Ml. 
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